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Semiconductor billiards are often considered as ideal systems for studying dynamical chaos in the quantum 
mechanical limit. In the traditional picture, once the electron's mean free path, as determined by the mobil- 
ity, becomes larger than the device, disorder is negligible and electron trajectories are shaped by specular 
reflection from the billiard walls alone. Experimental insight into the electron dynamics is normally ob- 
tained by magnetoconductance measurements. A number of recent experimental studies have shown these 
measurements to be largely independent of the billiards exact shape, and highly dependent on sample-to- 
sample variations in disorder. In this paper, we discuss these more recent findings within the full historical 
context of work on semiconductor billiards, and offer strong evidence that small-angle scattering at the sub- 
100 nm length-scale dominates transport in these devices, with important implications for the role these 
devices can play for experimental tests of ideas in quantum chaos. 
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1 Introduction 

Quantum mechanics and chaos theory are two great discoveries of the 20 th century. Combining them 
leads to highly interesting and important questions; for example, how does the uncertainty principle's 
inherent 'blurriness' affect the sensitivity to initial conditions characteristic of chaotic systems? Quantum 
chaos rose to popular attention in the 1970s, first in the context of electron dynamics in atoms 0]|2), 
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and subsequently for model two-dimensional (2D) dynamical systems known as 'billiards' |3][4). From a 
classical perspective, a billiard contains a particle traveling in straight lines between specular reflections at 
a shaped boundary wall. A variety of billiard shapes have been studied ranging from simple geometries 
(e.g., circle or square) to more complex geometries (e.g., stadia and iemon' billiards Q) and geometries 
with multiple topologically-independent walls (e.g., the Sinai billiard [§)). Billiards have since become an 
archetypal system for studies of quantum chaos QUO- 

The quest for experimental realizations of quantum billiards naturally followed, with microwave bil- 
liards formed in open resonators I19HTTI and superconducting slabs 1121 leading the way, closely followed 
by semiconductor microstructures lfP3}fi"5l . In the years that followed, quantum billiards have also been 
realized in quantum 'corrals' defined by scanning tunneling microscopy on metal surfaces |[T6l[T7l . in 
acoustic systems [18 19|, atom-optics systems |2() 21 1. and various optical systems including micro-cavity 
lasers 122] and double-clad optical fibers H23II241 . 

An important consideration in comparing measurements to theoretical expectations is to adequately 
account for practical aspects of the experimental implementation that affect the correspondence. For semi- 
conductor billiards, two important aspects are the electrostatic nature of confinement and impurities in 
the semiconductor 11251 . These aspects, impurity scattering in particular, were considered in initial stud- 
ies llT3lrT4ll2"6ll27l and suggested to not obscure the essential physics. This has resulted in a widely-held 
perception that semiconductor microstructures are an ideal test-bed for dynamical quantum chaos; for 
example, "... it is important to understand the consequences of non-diffusive electron dynamics on the 
electronic conductance or other transport properties. This question has been studied in much detail for 
semiconductor nano structures in which the motion of electrons is ballistic rather than diffusive. In such 
systems, disorder is negligible, and, consequently, all transport properties are determined by the shape of 
the sample, as in a billiard model." l28l . However, a range of experiments following these initial studies 
have gradually shown that impurity scattering not only cannot be ignored, but may actually dominate over 
boundary-dependent 'ballistic' transport in semiconductor billiards. 

Here, we provide a short, focussed review of studies of electron transport in semiconductor microstruc- 
tures, and offer evidence that these devices are not the ideal embodiment of quantum dynamical billiards, 
as is often thought based on the pioneering papers on this topic [13,14, 29 - 3 1 1 . Our aim is to show that the 
true picture of transport in these devices is much more complex than initially anticipated. 



2 Semiconductor quantum dots - The fundamentals 

2. 1 Making quantum dots 

The underpinning structure for most studies of semiconductor quantum dots and billiards is the modulation- 
doped AlGaAs/GaAs heterostructure 0321 . It typically consists of a GaAs substrate supporting an un- 
doped AlGaAs spacer layer, a Si-doped AlGaAs modulation-doping layer, and an undoped GaAs cap layer 
(Fig. 1(b)) grown using molecular beam epitaxy (MBE) 11331 . The narrow, triangular potential well at 
the AlGaAs/GaAs interface forms self-consistently through the band-bending induced by ionization of the 
Si dopants (Fig. 1(c)). The well width is comparable to the electron Fermi wavelength (Xp ~ 50 nm), 
strongly confining the electrons in the growth direction to form a two-dimensional electron gas (2DEG). 
Further confinement of the 2DEG is typically achieved using a set of metal surface electrodes (gates) de- 
fined using electron beam lithography techniques [34, 33 (Fig. 1(a)). Application of a negative bias V g 
to the gates electrostatically depletes the 2DEG directly underneath, transferring the gate pattern into the 
2DEG as shown in Fig. l(b)Q The result is a quasi-zero-dimensional quantum dot connected by small 
one-dimensional (ID) entrance and exit ports to 2DEG source and drain reservoirs. These are connected 
to the external circuit via diffused AuGeNi ohmic contacts. 



A common alternative is to selectively etch away the doping, which also gives patterned depletion regions in the 2DEG. 
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(a) • (b) (c) 

AIGaAs 
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Fig. 1 (a) Top view of a semiconductor billiard including basic measurement circuit. The billiard is defined by a 
pair of metal surface electrodes (gates) with a negative d.c. bias V g applied, and connected to the measurement circuit 
via 2DEG source and drain reservoirs, each featuring an AuGeNi ohmic contact. The dot-dashed line indicates the 
cross-section presented in (b). (b) heterostructure cross-section illustrating electrostatic depletion due to the negatively- 
biased gates, (c) Conduction band diagram illustrating 2DEG formation at the AlGaAs/GaAs interface. The dashed 
line indicates the electron Fermi energy Ef and + the ionized Si dopants in the n- AIGaAs layer. 



2.2 Length scales and transport regimes 

A hallmark of MBE-grown AlGaAs/GaAs heterostructures is an extremely high electrical mobility, typi- 
cally /i ~ 10 5 — 10 7 cm 2 /Vs, corresponding to an electron mean free path £ ~ 0.5 — 50 /im at a typical 
2DEG electron density n ~ 10 11 cm~ 2 . The traditional paradigm of transport regimes in mesoscopic 
devices ll36ll depends on the relationship between I and the device's length L and width W (Fig. 2(a-c)). In 
the diffusive regime, £ < L, W, the electron's trajectory is primarily determined by the impurity configura- 
tion (Fig. 2(a)). In the ballistic regime, £ > L, W, there are no impurities within the device on average, and 
an electron follows straight line paths between specular reflections at the device boundary (Fig. 2(c)). In 
the quasi-ballistic regime, £^L,WotW<£< L, a mixture of impurity and boundary scattering deter- 
mines the electron's trajectory (Fig. 2(b)). The quantum dots used as billiards typically have L, W < 1 /im 
and I > 2 /im. 

For later discussion, it is worth bearing in mind that this is a somewhat simplistic and incomplete picture 
of disorder in these devices, because /j, is heavily weighted towards large-angle scattering events 071 . As 
a result, small-angle scattering is often under-appreciated; we will return to this in §4.2. 





L ( 

Fig. 2 (a-c) The traditional paradigm of transport regimes in mesoscopic devices 1361 : (a) diffusive £ < L, W, (b) 
quasi-ballistic £ ~ L, W, and (c) ballistic £ > L, W. L and W are the device/channel length and width and £ is the 
electron mean free path determined from the electrical mobility /t. (d/e) Schematics for a modified picture of the (d) 
diffusive and (e) ballistic regimes accounting for small-angle ionized impurity scattering, as discussed in §4.2. The 
lighter dots in (d/e) correspond to small-angle scatterers, with darker shades representing a greater propensity towards 
larger scattering angles. 
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Fig. 3 (a) Scanning electron micrograph of a typical quantum dot device with Ti/Au gates (grey) on the GaAs surface 
(black). The lxl /im quantum dot is connected to 2DEG source and drain reservoirs via ID quantum point contact 
(QPC) entrance and exit ports, (b) Conductance G versus magnetic field B perpendicular to the 2DEG from the device 
in (a). Three traces are shown; the upper and middle traces are subsequent sweeps of B, the lower trace was obtained 
10 hours later, on returning the device to the same V g used for the first two traces after other measurements. 

2.3 Quantum interference and conductance fluctuations 

The Aharonov-Bohm (AB) effect (38l demonstrates the influence of electromagnetic potentials on quan- 
tum interference. The original concept involved an infinite length solenoid, with an electron beam split 
and recombined in the plane perpendicular to the solenoid axis to form an interferometer 'loop' . For an 
infinitely long solenoid, the field B experienced by the electrons is zero; the AB effect is caused by the 
magnetic vector potential A, such that electrons traversing opposite arms of the loop accumulate phase 
terms of opposing sign ±{e/h) J A ■ dl, where dl is an element of path length. The two electron beams 
interfere when recombined, giving a transmitted beam intensity that is an oscillatory function of the flux 
cj) enclosed by the loop (and in this case, restricted entirely to the solenoid's interior). Early electron mi- 
croscope experiments demonstrated flux-induced electron interference J39), but it took another 25 years to 
demonstrate interference occurring via the vector potential A alone [40 1 Q 

The Aharonov-Bohm effect lays at the heart of a range of magnetic field induced fluctuations in the 
electrical conductance of mesoscopic devices. At its simplest, it produces a sinusoidal fluctuation with a 
period AB = h/eA in the conductance G of nanoscale ring structures enclosing an area A; observed in 
both gold rings BTl and etched ring-shaped structures in an AlGaAs/GaAs heterostructure 11421 . These 
experiments are performed at a temperature T < 1 K using a 3 He cryostat or a 3 He/ 4 He dilution refriger- 
ator to ensure electron phase coherence. The conductance fluctuations develop richer structure when the 
interference involves a larger number of possible electron paths through the device. A classic example is 
the 1 /im square quantum dot in Fig. 3(a), which is connected to 2DEG source and drain reservoirs via a 
pair of Quantum Point Contact (QPC) entrance and exit ports. The QPC strongly diffracts the incoming 
electron wave [43 44], allowing a range of trajectories linking the entrance and exit QPCs to interfere and 
contribute to G. In a semiclassical picture, each trajectory pair encloses a different area 9, giving a slightly 
different period of AB oscillation; these combine to yield reproducible multi-spectral magnetoconductance 
G(B) fluctuations (Fig. 3(b)). These fluctuations are not noise; as the lower trace in Fig. 3(b) demonstrates, 
identical fluctuations are obtained if a device returned to a given gate bias configuration after measurements 
at other biases, providing the temperature remains low throughout (see §4.2). The link between the spec- 
tral content of the magnetoconductance fluctuations (MCF) and the distribution of loop areas N(9) R31 
provides a 'magnetofingerprint' of the electron trajectories within the dot 1461 1471 . While the underlying 
mechanism for these fluctuations is the same as that driving the well-known Universal Conductance Fluctu- 
ations (UCF) Il47l453l . the amplitude of the fluctuations for dots is not universal JT3); this is due to the lack 



Superconducting films were used to screen the electron beams from any magnetic field B generated whilst establishing A. 
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of ensemble averaging that occurs in diffusive systems [48. 49 52, 55] . Another commonly held difference 
between UCF in metals or MOSFETs and MCF in quantum dots is that whereas in the former transport is 
diffusive, in the latter the transport is assumed to be ballistic and hence determined by the geometry of the 
dot walls. We will return to examine this assumption in §4.2 in light of recent experiments [54. 55] that 
question the way we think about semiconductor billiards as systems for studying quantum chaos. 

3 Early pictures of quantum chaos in semiconductor billiards 

The use of ballistic microstructures for studies of quantum chaos 0561 evolved from experimental studies of 
ballistic junction devices 05711581 in the late 1980s. In particular, a seminal theoretical paper by Jalabert et 
al 11291 proposed that knowledge of the chaotic classical scattering dynamics could be used to predict quan- 
titatively measurable properties of ballistic microstructures. This relies on a universal property of chaotic 
dynamics in open scattering systems whereby a particle's probability of escape is exponential in time 059] , 
giving measurable correlations in the energy/frequency spectrum (60|6T|. Jalabert et al extended this prop- 
erty to ballistic microstructures, translating the escape time distribution into a distribution of accumulated 
trajectory areas, which, using semiclassical arguments, produces MCF with particular statistical properties 
via the AB effect. In particular, the autocorrelation function C(AB) =< 5g(B)Sg(B + AB) > for the 
MCF should take the form: 



where Sg(B) is the fluctuation from mean conductance. Here 4>q = h/e and a is the exponent for the 
distribution N(6) oc cxp(— 27ra|#|) of areas 6 enclosed by various possible trajectories 01311291 . 

Marcus et al lfT51 reported the first experiment devoted specifically to chaotic scattering in ballistic 
microstructures, measuring two separate chips, each containing two quantum dots, one with a circular ge- 
ometry and one with a stadium shaped geometry. The study had two intentions: to confirm the prediction 
by Jalabert et al 11291 for chaotic scattering systems, and to confirm predictions of quantitative statistical 
differences in the fluctuations between geometries supporting chaotic (stadium) and non-chaotic (circle) 
transport 06211631 . Comparison with theoretical expectations in Ref. |29l was made via the Fourier trans- 
form of Eq. 1, which gives the MCF power spectrum lfl"3l : 



Fits of Eq. 2 for the stadium geometry gave acf>o ~ 3.6 mT, corresponding to an a -1 ~ 2.7 times the dot 
area A, consistent with expectations from Ref. |29l . The spectral power S(f) for the circle was similar at 
/ < 250 T _1 but higher at 250 < / < 1000 T _1 , consistent with predictions of enhanced higher frequency 
content, corresponding to a larger number of large-area trajectories in non-chaotic geometries 0621631 . Note 
that the power spectra used for comparison to Ref. ||29l were heavily averaged (15 fast Fourier transform 
(FFT) spectra from half-overlapping 256-point (18.6 mT) wide windows). Single FFT spectra for 1024 
point blocks spanning —70 < B < 70 mT revealed structure at low frequencies, consistent with structure in 
the large AB tails of the C(AB) versus AB data presented in Ref. 03]. Such low frequency spectral peaks 
are an expected signature of wavefunction scarring by periodic orbits [4 44 64], and have been extensively 
studied in subsequent calculations 1 43 .44 .65 ]. transport studies 166) and scanning-gate microscopy (SGM) 
studies [67—69]. However, of particular note for the discussion in §4.2 is a comparison between the power 
spectra for nominally identical geometries (i.e., the stadium on device 1 with the stadium on device 2 or the 
circle on device 1 with the circle on device 2) in Fig. 3 of Ref. fT3l . If the dominant factor in determining 
the trajectory distribution is device geometry, then the spectra for identical geometries on different chips 
should match, or at least be very similar; they instead differ significantly, indicating that dot shape alone 
cannot account for the electron dynamics. 



C(AB) = C(0)/[1 + (AB/aM 2 } 2 



(1) 



S(/) = S(0)[l + (2Wo)/]e 



-27TQ0 O / 



(2) 
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Subsequent theoretical work by Baranger et al BUI proposed the alternative possibility of distinguish- 
ing between geometries supporting chaotic and regular (i.e., non-chaotic) dynamics using the 'ballistic' 
analog lfl3l of another important interference effect in diffusive 2D systems known as weak localization 
(WL) l70l . Weak localization involves interference between a trajectory that returns to its point of origin 
via a set of impurity scattering events, and its time-reversed counterpart. Application of an increasing per- 
pendicular magnetic field B progressively breaks time-reversal symmetry for shorter 1 ength paths, leading 
to a positive magnetoconductivity correction (minima in G(B) at B = 0), and hence a peak in the mag- 
netoresistance.[] The proposal by Baranger et al 1301 is that the different dynamics produce different WL 
peak lineshapes, which relies on the distribution of enclosed loop areas N{&) differing between chaotic 
and regular billiards. Under a semiclassical framework, the average quantum correction to the reflection 
coefficient Rp for chaotic billiards takes the form: 



for small B BUI . It is important that the average above is taken over an infinite window in k to ensure that 
the contributions to Eq. 3/4 only include the interference of paths that are precise time-reversed pairs. This 
averaging has important experimental implications - the geometry effect should not be directly observable 
in a single magnetoresistance R(B) trace, but only when sufficient averaging is applied to eliminate the 
conductance fluctuations l74l . Indeed, it is reported in Ref. l74l that calculations of unaveraged R(B) 
traces yield a minimum at B = in ~ 33% of cases, in contrast to the expected B = maxima. This is 
borne out experimentally; e.g., in studies of dots featuring an additional 'shape-distorting' gate by Chan et 
al 1751 averaging was needed to obtain a clear Lorentzian R(B) peak and R(B) minima are occasionally 
observed (76], in studies of a square billiard by Bird et al 1771 where R(B) peak lineshape (Lorentzian 
vs linear) was sensitive to gate voltage, and in studies by Keller et al 11781 where G(B) maxima and 
minima were obtained at different densities (i.e., Fermi energy) while the ensemble-average (G(B)) gave 
the expected minima at B = 0. 

Measurements by Chang et al lfT4l provided evidence in support of the prediction by Baranger et al 
above; reporting a Lorentzian R(B = 0) peak for a stadium geometry (chaotic) and a linear R(B = 0) peak 
for a circular geometry (non-chaotic). It is very important to note that these studies were not performed 
in single dots, but in devices consisting of 48 dots with nominally-identical geometry arranged in a 6 x 8 
array. This is vital in terms of the averaging issues discussed in Ref. l74l . but may also have important 
implications in terms of small-angle scattering, as we will discuss in §5. 

This brings us to an important turning point in terms of our understanding of ballistic transport and 
quantum chaos in semiconductor billiards and quantum dots. On the basis of the work discussed above it 
is tantalizing and easy to assume a great simplicity in these devices: images of electrons traveling along 
straight line trajectories between specular reflections from walls with well defined geometrical shapes, 
clear divisions of geometries into chaotic and non-chaotic, which in turn give clear and distinct experi- 
mental signatures, etc. see, e.g., Ref. PP . As the following discussion will show, in reality, the picture of 
electron dynamics in semiconductor billiards is more complex. The practicalities involved in the experi- 
mental realizations have effects that are far from negligible, and that will ultimately force us to rethink our 
understanding of the results previously described in this section. 



n.b. alternative explanations for the magnetoconductivity about B = in quantum dots have also been subsequently pro- 



(5R D {B))=R/[l + {2B/a<t> ) 2 \ 



(3) 



where R = (5R D (B = 0)} El, and for regular billiards: 



(SR D (B)) oc \B\ 



(4) 



posed 17 141731 . 
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Fig. 4 (a) Gate schematics (top) and scanning electron micrographs (bottom) for the seven billiard devices studied in 
Ref. 1851. Devices E, F and G are identical to B with with widths of 1, 0.6 and 0.4 /im, respectively. The digits in (a) 
indicate lengths/widths in [im\ Ti/Au gates appear grey on the GaAs surface (black), (b) Measured fractal dimension 
Df versus Q, which quantifies the resolution of quantized energy levels within the billiard 1851 . (c) Fractal scaling 
plot for G{B) data from the 'disrupted' billiard from Ref. J89] at T = 50 mK (top), 500 mK and 1.2 K (bottom). The 
arrows indicate the lower cutoff for fractal scaling. 



4 The trouble with semiconductor billiards - Soft walls and small-angle 
scattering 

4.1 Soft-walls and mixed phase-spaces 

The complexities involved in the theoretical calculations described above typically forced semiconductor 
billiards to be stripped down to their most basic elements for modeling (particularly in the early 90s, where 
computational capabilities were not what they are today). The early models of dynamics in semiconductor 
microstructures assumed hard-walls (i.e., the location of the wall in terms of particle reflection is indepen- 
dent of energy) and geometries where the phase space is either completely chaotic (ergodic) or completely 
regular. However, semiconductor billiards are defined by electrostatic depletion of selected regions of the 
2DEG, using either negatively biased surface-gates or by locally etching away the dopants. This naturally 
results in a soft-wall potential 17911801 . often best approximated as parabolic (25). This profile results 
generically in a mixed phase-space, containing both chaotic and regular regions; the ramifications of this 
were explored by Ketzmerick l8D . Trapping by the infinite hierarchy of Cantori at the boundaries between 
regular and chaotic regions in the phase space causes the escape time probability distribution to take a 
power-law form P(t) ~ t~P with exponent f3 < 2. This was expected via semiclassical arguments to pro- 
duce a power-law distribution of enclosed areas N(6) ~ 9 1 for large 9 and statistically self-similar fractal 
conductance fluctuations with fractal dimension Dp = 2 — 7/2 for 7 < 2 (and Dp = 1 for 7 > 2) (STJ. 

The observation of statistically self-similar fractal conductance fluctuations in billiards was indepen- 
dently reported first by Micolich et al l82l and soon thereafter by Sachrajda et al (83)Q In both cases, 
fractal dimensions 1 < Dp < 2 were reported, with preliminary analyses reporting dependencies of Dp 
on external parameters such as gate voltage and temperature 18211831 . 

A subsequent, more extensive study of the dependence of Dp on a range of parameters involving seven 
separate devices (Fig. 4(a)) by Micolich et al (85) lead to a remarkable finding: a plot of Dp versus a 
single parameter Q for many G(B) traces all fall onto a single curve (Fig. 4(b)), despite being obtained 
from seven different sizes and shapes of device with differing number of ID modes n passing through the 



n.b., Exact self-similarity in a semiconductor Sinai billiard was reported by Taylor et al 1 84 ] a year prior. 
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QPC entrance and exit ports, electron mean free path I and temperature T. The parameter Q is the ratio of 
the average energy level spacing AE$ to the average energy broadening AEb for the dot, taking the form: 

AE S _2^ 1 

AE B m*A ^{h/ Tq Y + (k B T)^ 

where m* is the electron effective mass and r q is the lifetime of the quantum states, obtained using a skip- 
ping orbit analysis 11861 488 1 1^ Remarkably, Dp is found to take its maximum value ~ 1.52 at Q = 1, with 
Dp falling in either direction, reaching the non-fractal value Dp = 1 in the limits Q = and Q ~ 10. 
The trend in the small Q limit is particularly interesting; small Q is typically achieved via the denomi- 
nator through increases in T and/or decreases in r q . Based purely on Ref. JSTJ, one would expect that 
a reduction in the coherence length (and thus r q ) independent of other billiard parameters would reduce 
the scaling range over which fractal behaviour was observed but not alter Dp itself. Indeed, later work 
(Fig. 4(c)) showed that the scaling range instead increases with increasing temperature concurrent with a 
falling Dp ||89l . Calculations by Budiyono and Nakamura [90 1 based on a semiclassical Kubo formalism 
obtained a Dp that decreases smoothly with increasing T; an interesting side-aspect of these calculations 
is a connection to the Weierstrass function model [91. 92) for the exact self-similarity reported in Sinai 
billiards 11841 . Hennig et al 11931 provide an interesting alternative, suggesting that the underlying classical 
dynamics also lead to fractal conductance fluctuations. As such, G(B) is a superposition of two fractals, 
one due to quantum interference, which is suppressed by decoherence to reveal fractal fluctuations due 
to the classical processes. Note that the data below the lower cut-off (i.e., low AB limit) in Fig. 4(c) is 
linear with a different slope. This might be indicative of a process similar to that proposed by Hennig 
et al, however, care is needed as this data is close to the resolution limit/noise floor; further experiments 
in this direction are warranted. Several studies suggest that the semiclassical theory proposed by Ketzm- 
erick (STj is only part of a larger picture. A complete discussion is beyond the scope of this paper, but 
most significantly, mixed phase space does not appear essential to the generation of fractal conductance 
fluctuations [94 , 95 1, which are also expected theoretically for completely chaotic [95 -97 1 and completely 
integrable (98) billiards. 

Considering Ketzmerick's theory, an interesting parameter to vary is the softness of the billiard walls. 
The first studies were performed by using a pair of billiards defined, one above the other, in a double quan- 
tum well heterostructure to provide identical geometries with different wall softness 11991 . Self-consistent 
potential simulations 11 001 showed that the profile gradient at the Fermi energy differed by a factor of 
three between the 90 and 140 nm deep billiards 1991 . The data in both billiards was consistent with Eq. 5; 
the only difference compared to the data in Fig. 4(b), which contains the shallower billiard data as device 
C, is that the deeper billiard (harder potential) gave a lower maximum Dp <~ 1.45 at Q = 1 l89l . This 
difference in maximum Dp is small, and may be due to competing effects given the complexities of this 
device l99l . 

The effect of soft-wall profile was also studied by Marlow et al 11011 . A further 21 devices were 
investigated: 11 Gao.25Ino.75As/InP billiards and 5 Gao.25Ino.75As/InAlAs billiards of various geometries, 
and 5 GaAs quantum wires. As demonstrated by Martin et al 11021 . the potential profile for the etched 
GalnAs/InP billiard^ is ten times steeper than for the surface gated AlGaAs/GaAs billiards in Ref. 1851891 . 
Despite this, the data for all of these devices closely overlays that in Fig. 4(b), suggesting that wall profile 
is not a determining factor for Dp. The data obtained for quasi-ballistic and diffusive GaAs quantum 
wires was also fractal and fit on the curve in Fig. 4(b). This is consistent with earlier work showing fractal 
conductance fluctuations in quasi-ballistic gold nanowires 11031 . and predictions that fractal conductance 
fluctuations survive diffusive transport 11971 . 



This is ultimately an unentangleable combination of the phase coherence time and the escape time related to a skipping orbit 
finding a QPC through which to exit the dot. 

^ The substitution of InAlAs for InP should not affect the profile of the billiard walls. The slight change in band-gap (from 
~ 1.42 eV to ~ 1.55 eV) should not significantly affect the 2DEG confinement either. There is a large difference in disorder 
however, as described in the text and Ref. 11011 . 
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Ultimately, the results obtained by Micolich et al 118511891 and Marlow et al 11 1 II suggest that the 
statistical properties of the conductance fluctuations, at least in terms of fractal dimension, are independent 
of the geometry of the billiard wallfl This is unexpected given that transport in these devices is nominally 
'ballistic'. Instead, the statistical properties via Eq. 5 and Q, depend only on A, T and r q irrespective 
of the specifics of the sample itself. This 'universality' bears a strong resemblance to that of Universal 
Conductance Fluctuations (UCF) in mesoscopic metal films and MOSFETs (see, for example, Fig. 3 of 
Skocpol 153]) where the transport is diffusive rather than ballistic. This raises interesting and radical 
questions: Even though I » L, W for billiards and they are ballistic by traditional measures 1361 . is the 
transport really ballistic? Or is it diffusive, such that a form of scattering essentially undetected by the 
mobility, dominates the transport? 

4.2 Remote ionized dopants and small-angle scattering 

The 2DEG in an AlGaAs/GaAs heterostructure is normally populated by the ionization of remote Si 
dopants. These positively charged ionized dopants interact electrostatically with the electrons in the 2DEG 
causing significant scattering. High mobility 2DEGs are obtained through a process known as 'modulation 
doping' l32l . where the Si dopants are separated from the 2DEG by a 20 — 100 nm undoped AlGaAs spacer 
layer. The key to modulation doping is that ionized dopants are converted from large-angle scattering sites 
to small-angle scattering sites by increases in the dopant-2DEG separation. This effectively 'hides' some 
of the ionized dopant scattering because the mobility is weighted towards large-angle scattering D7II1051 
(see §5). Nonetheless, the 2DEG still feels the effect of the ionized dopants, which present as a low-level 
(~ 1 — 10 meV) random 'disorder potential' for the electrons [106]. The length scale of this scattering is 
set by the 2DEG-donor separation 111071 . usually of order 20 — 100 nm, is much smaller than the typical 
billiard width (~ 0.6 — 2 ^m) and the large-angle scattering length £ ~ 2 — 20 /im. 

At the time of the initial studies II 13 II 141 29 30 36 1 , knowledge of the effect that the underlying disorder 
potential had on ballistic transport in microstructures was relatively low, although numerical simulations 
showed it could significantly affect transport at length scales much smaller than I | 108 1 . The first focussed 
consideration in terms of billiards was by Lin, Delos and Jensen 11271 . and involved adding a random 
angle [— tv, +tt] at each wall reflection in a semiclassical dynamical model for stadium and circle billiards. 
The focus was on the statistics, and while the disorder in this model had little effect for the stadium (as 
expected given the ergodic dynamics), it destroyed the power-law tail in the area distribution for the circular 
billiard. This gave an exponential distribution, as expected for the stadium billiard, albeit with a different 
exponent 11271 . This should produce a Lorentzian R(B = 0) peak for both stadium and circle billiards. 
Chang et al lfl4l approached the same problem differently, adding a random energy [— Wdis/2, Wdis/2] 
at every fifth lattice site and linearly interpolating in between in quantum calculations on a discretized 
lattice 01091 . They concluded that the linear R(B = 0) peak for the circular billiard should not be destroyed 
by a disorder potential with a strength chosen to match experimental conditions, with roughness in the dot 
wall being more crucial lfT4l . Ji and Berggren Il25lll 101 made the first investigation at the single G(B) 
trace level (as opposed to the statistics of area distributions or averaged R(B = 0) peak lineshape) for 
stadium-shaped billiards using a discrete lattice model 111 1 II . Their model showed that changes in disorder 
potential can have a significant impact on G(B) as a magnetofingerprint for the electron dynamics [1 10 1, 
as predicted for UCF by Feng, Lee and Stone 11461 . and demonstrated for diffusive quantum wires by Taylor 
et al 111 121 and Klepper et al 11 131 . 

Recently, a technique known as scanning gate microscopy (SGM) Ml 141 has given remarkable insight 
into the 'disorder potential' and its influence on electron transport 11 1 0711 1 15lll 161 . SGM studies of 2DEGs 
show significant branching of electron flow 111 151 . and even in samples with very long mean free paths 
(£ > 10 /im) significant deviations in propagation persist at scales ~ 100 nm 110711 . An SGM study of a 
ballistic focussing structure by Aidala et al 11 161 gives a remarkable visual demonstration of the impact that 
small-angle scattering has on even quite short and simple ballistic transport paths between QPCs. Although 



Note that other statistical measures of the fluctuations, such as the power spectrum, are also linked to fractal dimension 11041 . 
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small-angle scattering clearly doesn't prevent ballistic focussing from being observed electrically 01 171 . it 
certainly forces one to question whether the concept of semiconductor billiards with long straight-line 
trajectories bouncing specularly off the billiard walls (see e.g., Fig. 2 of Ref. OTI or Fig. 1/6 of Ref. [5J) 
is a realistic picture? It also raises the question of the extent to which small-angle scattering governs the 
dynamics and resulting electrical properties of these devices. 

Scannell et al |54j and See et al 11331 have recently performed experiments aimed at illuminating the 
relative importance of disorder and boundary scattering in semiconductor billiards. These experiments 
rely on an important feature of the dopant physics of AlGaAs/GaAs heterostructures and the newly de- 
veloped capability to make billiards where the 2DEG is populated electrostatically using a gate enabling 
the dopants to be removed 111 181 and the effect of small-angle dopant scattering on transport to be clearly 
observed. Si dopants in modulation-doped Al^Gai-^As/GaAs heterostructures with x > 0.22 can take 
one of two possible configurations 01 19111201 . The first is as a shallow hydrogenic donor with a positive 
or neutral charge (d + /d°) sitting at a substitutional site. The second, known as a DX center, involves a 
lattice distortion that pulls the Si donor into an interstitial position along the (111) direction, producing 
a deep trap. At temperatures below ~ 150 K, the DX centers normally bind two electrons and have a 
negative charge (DX~), but can be excited thermally or optically to detrap into DX° and DX + states. As 
described earlier, the ionized dopants impose a 'disorder potential' on the 2DEG 1 106|, with a density of 
scattering sites N = N(d + ) + N(DX ) + N(DX + ), where N(...) is the density of a particular dopant 
configuration. While N is normally fixed for a given heterostructure and gate-voltage configuration during 
cooldown II 1201 . the allocation of particular Si site as d or DX varies varies each cooldown. Thus the result- 
ing disorder potential is fixed for a single cooldown, leading to the repeatability of G(B) traces shown in 
Fig. 3(b) and Ref. 11261 . In contrast, the disorder potential should change markedly on separate cooldowns, 
and this would produce one of two results. If the transport in the billiard is truly ballistic, such that the 
disorder is negligible and the electron trajectories are determined only by scattering from the boundaries, 
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Fig. 5 Conductance G versus magnetic field B and gate voltage V from (a/b) a modulation-doped billiard and (c/d) 
an undoped billiard. False colour scanning electron micrographs of these devices are shown inset to (b) and (d) 
respectively. The scale bars indicate 500 nm. The two devices have I = 2.7 and 2.1 ^im respectively. In (b), the gate 
bias V g is applied to both metal gates (yellow). In (d), the gate bias Vtg is applied to the n + -GaAs top-gate (green), 
with the side-gates (red) grounded. In both cases B = 0. The traces in (a) were obtained at V g = —430 (top) and 
—451 mV (bottom); these voltages were chosen to match G(B = 0) before and after cycling. G(B) traces obtained 
with common V g do not match either (see Ref. 1551 ). The three pairs of traces in (c) were obtained at Vtg = +955 
(top), +945 (middle) and +935 mV (bottom). The traces are shown in each case from top down: first (black), second 
(red) and third (blue) cooldown. Traces are vertically offset for clarity by: (a) bottom trace —0.3; (b) top trace +0.9, 
middle trace +0.3; (c) from bottom by -0.3, 0, +0.7, +1.0, +1.7, +2.0; (d) top trace +0.3; all in units of 2e 2 /h. 
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then G(B) should remain the same in subsequent cooldowns (providing T > 150 K in between). Alterna- 
tively, if the disorder potential does affect transport, then we should see some level of change in G(B), as 
predicted by Feng, Lee and Stone lt46l . 

Figure 5(a) shows G(B) traces obtained on separate cooldowns of the modulation-doped device inset 
to Fig. 5(b). There is clearly a striking difference between the traces indicative of the disorder potential 
affecting the transport. The gate voltage was adjusted slightly on the second cooldown in order to match 
G(B = 0) and obtain similar injection conditions 043 11441 . but See et al observe the same behaviour if V g 
is kept the same, or takes other values as part of an extensive search for matching G(B) traces (55] . As 
Fig. 5(b) shows, the change in disorder potential upon thermal cycling also alters the gate characteristics; 
the fluctuations there also originate from quantum interference (66j . This thermally-induced irreproducibil- 
ity is also evident in a parallel study by Scannell et al 1541 . where the device is cooled to 300 mK and G(B) 
is measured, the temperature is raised to an intermediate temperature T\ for 30 minutes and returned to 
300 mK, G(B) is measured again, the temperature is raised to T-2 > Xi, re-cooled and measured again, 
and so on until T n = 298 K. The various post-cycling G(B) traces are compared to the first cooldown 
G(B) trace using a cross-correlation analysis 111211 . with a precipitous drop in the correlation coefficient 
F from 1 to obtained at T — 100 - 175 K for AlGaAs/GaAs and GalnAs/InP billiards and GaAs quantum 
wires. This temperature dependence points strongly to ionized dopant scattering playing a major role in 
influencing transport in these devices. 

The 'acid test' for the hypothesis that the disorder potential dominates transport is to remove the dopants 
entirely. This requires a complete redesign of the device; a metallic gate needs to be added and this needs 
to overlap the contacts to the 2DEG to ensure electrical continuity. This is achieved using a degenerately- 
doped GaAs cap layer as the gate, which enables self-alignment of the gate to the contacts 11 18111221 . The 
n + -GaAs cap is divided into three gates by shallow wet etching as shown inset to Fig. 5(d); the central gate 
(green) is biased positively to Vtg to electrostatically populate the 2DEG, the side-gates (red) can be used 
to tune the device further, but remain grounded here. Figure 5(c) shows G(B) before and after thermal 
cycling to 300 K for three different Vtg- m eacn case the MCF are reproducible with high fidelity. As 
Fig. 5(d) shows, the gate characteristics are also highly reproducible. 

Given the MCF's role as a 'magnetofingerprint' lt46l of electron transport within the dot, Fig. 5 demon- 
strates that the disorder potential has a significant, perhaps even dominant, influence on the electron dy- 
namics. We now discuss an explanation for how and why this occurs, the implications regarding past 
results, and some thoughts on follow-up experiments to further elucidate the role that disorder plays in 
billiards. 



5 Discussion 

Historically, a major focus in mesoscopic physics has been to improve the mobility to access new physical 
regimes that are otherwise obscured by disorder. This was a key driver for work on the fractional quan- 
tum Hall effect 1123111241 and ballistic transport devices (e.g., quantum wires, dots and billiards.) as per 
Figs. 2(a-c). It is vital to note that the mobility, and by association I, are only part of the picture. The 
mobility is a measure of momentum decay and it contains a weighting term 1 — cos 9 1371 . where 8 is the 
scattering angle, that makes it preferentially sensitive to large-angle scattering events [37 105 107] Fig- 
ures 2(d/e) present a more refined picture for the diffusive and ballistic regimes; here large-angle scattering 
sites are black and small-angle scattering sites are grey (darker shades represent a propensity towards larger 
scattering angles). The two pictures of diffusive transport (a/d) are essentially the same, and in either case, 
it is impurity scattering that dominates over boundary scattering. The two pictures for ballistic transport 
(b/e) are very different; as in §4.2 the ionized dopants present as small-angle scatterers with a scattering 
length set by the spacer thickness [ 107 1, which is much smaller than L, W and I. These deflect the electron 
trajectories between the walls as shown, suggesting a more complex picture with a less clear-cut distinc- 
tion between the traditional diffusive and ballistic regimes. Bearing in mind the thermal redistribution in 
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dopant potential [54, 55 1, this explains why the MCF changes so much between cool-downs and differs 
when studying the same geometry on two separate chips (see Fig. 3 of Ref. Ifl3l ). 

From this perspective, the MCF in modulation-doped billiards is probably best viewed as analogous to 
UCF in metal films and narrow MOSFETs 047H531 . In both cases the transport is dominated by disorder 
scattering, the fluctuations are reproducible in a single sample/disorder configuration but vary between 
samples and disorder configurations, are quenched by an increase in temperature, and have certain univer- 
sal statistical properties independent of the sample shape and precise dopant configuration [47i| - for the 
UCF this is the amplitude 5G which depends on the ratio of length L to coherence length 05 31 . whilst 
for MCF in billiards, this appears to be the fractal dimension Dp which depends on a more complex com- 
bination of area A and quantum lifetime r q , at fixed temperature 0851 . The billiard walls should help rather 
than hinder this analogy; their action is to repeatedly feed the electron back into the diffusive environment 
inside the billiard 1 1011 . In a sense, the walls impose reflective periodic boundary conditions, effectively 
converting the billiard into a larger area, diffusive 2D system. The natural correlations between these mir- 
rored sub-units, combined with the two QPCs acting as points of focus for the interference contribution 
may explain the loss of amplitude universality, and its replacement by the more complex universality in 
Fig. 4(b) [85 89 101 1. Modeling a small-area, small-angle diffusive system with reflective boundary con- 
ditions on the edge for all but two points where trajectories can enter/exit and contribute to G(B), and 
comparing with traditional UCF 0471 may give interesting insights. Alternatively, it may be interesting to 
consider existing UCF models, but with the uniform random disorder replaced by repeating small 'tiles' of 
a common disorder configuration. 

Nonetheless, 'ballistic' transport, where the walls rather than disorder is the shaping force for a given 
trajectory, cannot be neglected entirely in billiards. Comparing Figs. 2(b,c,e) it is clear that the diffusive 
effect of small-angle scattering should increase with trajectory length. As a result, the long trajectories 
crucial to discriminating between chaotic and regular geometries in semiclassical theories O29"ll30l are 
likely destroyed by disorder such that all billiards are essentially chaotic, irrespective of geometry, from 
the perspective of the tails in the trajectory area distributions N{6). As such, the disorder has the same 
'chaos-inducing' effect as a Sinai diffuser l6l l6TIIT00l . but with a scale and quantity more similar to antidot 
arrays B 12511 1261 where chaotic transport, weak localization and quantum interference fluctuations are also 
observed. Shorter trajectories may well survive, and we see two separate categories of these. The first are 
skipping orbits as highlighted by Christensson et al 111271 . and used for studying phase-breaking in dots 
by Bird et al |f86l l88l . These are likely protected from disorder scattering by the short distance between 
bounces combined with a field-induced curvature mechanism similar to that involved in the quantum Hall 
effect 01281 . The feasibility of the survival of skipping orbits amidst the disorder potential is evident 
in the SGM studies by Aidala et al 01 161 . The survival of short ballistic trajectories may also explain 
measurements of the symmetry of MCF in the non-linear voltage regime as a function of the symmetry of 
billiard shape [129 1 30 1 . Combined, the work above provides strong evidence that billiard shape is still 
exerts an influence on the MCF. 

The second category is periodic orbits, the classic example being the diamond scar in square dots 0431 
[44); this also clearly survives small-angle scattering as observed in SGM studies by Burke et al 0681 . 
but scars with other shapes have also been reported very recently J69). It is easy to consider this as 
a straightforward case of wavefunction scarring [4. 44. 64 . but there may be more to it. The diffusive 
effect of disorder should not only affect the unstable periodic orbit underlying the wavefunction scar, but 
it will do so differently between devices, between different cooldowns of the same device, and possibly 
also as parameters such as the gate voltage are changed. The latter may partially explain the observation 
of some but not all of the expected scar states in previous experiments on modulation-doped dots 0661 . 
Wavefunction scarring in dots has recently been linked to a theory known as 'Quantum Darwinism' 01311 . 
where preferred states within the dot survive coupling with the environment, which acts as the 'selection 
pressure' in the theory 111 321 . From an experimental perspective 0661 . it may be that disorder presents 
an additional selection pressure, acting in partnership with environmental coupling, that eliminates scar 
states connected to much more complex underlying periodic orbits. From this perspective, it would be 
interesting to extend recent SGM studies of dots [68 [69] to undoped heterostructure system to see if a 
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greater variety of scars are observed. It would also be interesting to perform SGM studies of undoped 
2DEGs and ballistic focussing structures more generally to establish the level of improvement provided by 
eliminating the ionized dopants. The current architecture, with a gate above the 2DEG, is unsuitable for 
SGM; an alternative would be to use an inverted structure 11133111341 . possibly using deep etching instead 
of surface gates and/or patterning of the doped gate layer using focussed ion-beam lithography, as in the 
bilayer billiard ||99l . 

We conclude by attempting to 'close the loop' and reconsider ballistic microstructures as devices for 
studying dynamical chaos in the quantum limit. The data in Fig. 5(c/d) shows that removal of the modula- 
tion doping instils thermal robustness in the MCF obtained from the device. Correspondingly, one would 
expect a significant reduction in small-angle scattering, despite this, substantial disorder will remain. This 
includes background impurities, some ionizable, some not, as well as interface roughness, defects in the 
walls of the device, etc. We suspect that this remnant disorder will prevent perfectly identical G(B) from 
being obtained in separate undoped devices with nominally identical geometry, obviating the truly ballistic 
quantum dots envisioned theoretically 129H3TI . One path forward might be the approach used by Chang 
et al lfl4ll . where a 6 x 8 array of nominally identical dots was studied to average out the MCF. This av- 
eraging of the MCF essentially constitutes an averaging of the dot disorder potentials to expose the effect 
of the common lithographic geometry. Attempting this using arrays of undoped dots in clean heterostruc- 
tures with low background impurity densities may bring these devices close to a true 'ballistic' limit in 
terms of semiclassical dynamics, providing of course that the lithography defining their walls can be made 
sufficiently identical. This would be a useful challenge for future work. 

Ultimately, while single quantum dots are not sufficiently 'ballistic' for use in studies of quantum chaos 
in our opinion, they have been shown to be sufficiently ballistic for potential applications II 1 351 as ballistic 
rectifiers [136] and Y-branch switches [11371 . The reduced small-angle scattering in undoped quantum dots 
offers much potential in this regard, their thermal robustness suggests they may also offer a path to devices 
with significantly reduced charge noise, which is a well known effect of fluctuations in the charge state of 
ionized dopants 111381 . 
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